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Abstract 


The  phase-lock  phenomenon  of  a  mass  sled  sliding  along  in  a  circular 
slingatron  is  studied  both  numerically  and  analytically.  Parameters  that  describe 
a  slingatron,  in  which  the  phase  angle  of  die  swing  arms  increases  quadratically 
in  time,  are  found  to  be  simply  related  to  the  sled's  speed  during  phase  lock.  The 
time  required  for  phase  lock  to  occur  is  related  to  a  simple  exponential  function 
of  the  gyration  speed  and  the  coefficient  of  friction  between  tihe  sled  and  track. 
Accurate  time  histories  describing  the  motion  of  the  accelerating  sled  are 
expressed  in  terms  of  confluent  hypergeometric  functions  ,F, .  These  results  are 
then  used  to  obtain  physical  insight  into  why  the  phase-lock  phenomenon  takes 
place  and  to  describe  the  important  role  that  friction  plays  by  damping  the 
oscillatory  motion  of  the  sled. 
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1.  Introduction 


A  concept  called  the  shngatron,  used  to  accelerate  a  mass  (sled)  to  high  velocity, 
has  been  proposed  and  studied  by  Tidman  [1].  This  accelerator  can  propel  either 
a  constant  mass  or  an  ablating  mass  sled  along  a  guide  tube  forming  the 
slingatron  track  [2];  several  slingatron  configurations  have  been  examined 
[1, 2].  A  simple  friction  model  is  also  used  in  which  the  friction  force  is  assumed 
to  be  proportional  to  the  normal  force  exerted  by  the  track  [2]  on  the  moving  sled. 

Previous  work  has  shown  that  the  relative  phase  angle  (between  the  sleds' 
velocity  vector  and  the  gyration  velocity  vector  of  the  guide  tube)  locks  into  an 
approximate  constant  value.  This  report  presents  numerical  and  analytical 
results  foxmd  from  our  investigation  of  this  phase-lock  phenomenon  for  a 
constant  point  mass  accelerated  aroimd  a  circular  track  with  the  friction 
coefficient  taken  to  be  a  constant.  Numerical  simulations,  as  well  as  simple 
analytical  results,  show  that  the  mass  can  reach  very  high  velocities,  provided 
that  the  gyration  speed  can  increase  with  sufficient  acceleration.  We  will  show 
that  when  the  gyration  phase  angle  varies  quadratically  with  time,  tihe 
phase-lock  value  plus  the  time  required  to  obtain  this  value  are  simply  related  to 
the  gyration  acceleration,  the  coefficient  of  friction,  and  the  ratio  of  the  two 
characteristic  lengths  that  describes  the  slingatron  geometry.  These  relationships 
are  used  to  aid  our  understanding  and  to  make  predictions  of  the  phase-lock 
phenomenon. 


2.  Equations  of  Motion 


Consider  the  circular  slingatron  shown  in  Figure  1,  which  has  radius  vector  R  at 
the  sled  location  with  its  center  attached  to  the  gyration  arm  r  positioned  at 
angle  ^(t)  with  respect  to  the  laboratory  frame  at  time  t .  The  position  of  the 
accelerating  point  mass  is  given  by  R  +  r,  and  the  corresponding  velocity  is 
therefore  V  =  R  +  f .  The  force  F  acting  on  the  mass  M  can  be  written  as 
F  =  Fj^m  +  F//I1 ,  in  which  m  =  n  x  k  and  n  =  -  R'/R'  are  unit  vectors  pointing 
normal  (toward  the  center  of  the  circle)  and  anti-parallel  to  the  circular  track. 
Letting  i,  j  be  ;mit  vectors  along  the  laboratory  frame's  abscissa  and  ordinate 
gives  the  time  derivative  of  the  momentum  equations  the  following  form: 

Mx  =  -Fj^coscj)  +  F//sm(l)  and  My  =  -Fj^sincf)  -  F// cos<|) ,  (1) 
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Circular  Slingatron  Track 


Figure  1.  Schematic  of  a  circular  slingatron. 

with  Cartesian  coordinates  (x,  y)  in  the  (i,j)  laboratory  frame  constrained  sc 
that 

X  =  Rcos<|)  +  rcosvj/  and  y  =  Rsm<|>  +  rsinvj/ .  (2] 

The  parallel  the  force,  V,/ ,  attributed  to  friction  is  represented  cis 

Fii=t^P'x  •  (3) 

Equations  (1-3)  will  determine  the  angle  (|)  in  terms  of  the  given  angle  xi;(t) 
resulting  from  the  following  differential  equation 


-(^  =  +a[\j/  +  iu\t/^)cos{\j/- ^)+a[ny/-\ir^) sin  {y/-<l))  (4) 


for  which  the  length  ratio  is  defined  as  a  =  r/R . 

Generally,  tihe  solution  to  equation  (4)  is  obtained  from  numerical  integration  for 
any  choice  of  the  time-dependent  gyration  angle  X(;(t).  The  mass  velocity 
+  y^  relative  to  the  laboratory  reference  frame  is  easily  calculated  using 


V^  = 


2a  cos(v|/  -  (|))—  +  + 


xj;'  R' 


(5) 


3.  Phase  Lock  In 


Tidman  [1]  has  presented  a  complete  analysis  of  cases  where  phase  lock  in  occurs 
for  aU  times,  t ,  when  the  angle  0  =  v|/  —  (j)  is  constant,  so  that  ^  =  vj;  and  ^  =  xj/ . 
These  conditiorxs  transform  equation  (4)  to 


y  ,  ,  aism.0  -  /ucosO)-  jJ. 

(j)  — “  b  (p  ^  b  y  \  / 

ayjusmO  -¥cos6  J  +  1 


(6) 


and  this  is  easily  integrated  so  that  equation  (5)  becomes 

^2  _  (<^0  R)^(2aco5  0  +  a^  +i) 

with  <t>o  =  <t>  {o)  and  ({>0  =  <j)  (d). 


Examining  the  factor  b ,  defined  in  the  second  equation  of  equations  (6),  shows 
that  the  velocity  is  time  independent  when  b  =  0  so  the  parameters  are  related  as 


0  =  ±2  tan 


-1 


■yj (a^  +a^  ±a 

(a-7)p 


±  n  TT,  n  =  0,7,2  •  •  • , 


so  that  (8) 

Pl^+1 


A  short  amoxmt  of  algebra  reveals  that  equation  (8)  is  satisfied  whenever  the 
force  caused  by  friction,  p  Fj^ ,  is  large  enough  to  prevent  the  point  mass  from 
accelerating  in  a  direction  tangent  to  the  circular  track.  Furthermore,  the  velocity 
given  by  equation  (7)  was  shown  by  Tidman  [1]  to  diverge  (to  a  nonrelativistic 
limit)  at  the  time  t„  =  ^  practice,  for  tiiis  to  occur,  one  must  have 

t„  >  0,  which  means  that  a  divergent  velocity  V  wiU  occur  xmder  the 
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(9) 


assumption  a^Jl  +  ju^  <  1,  whenever  the  parameters  also  satisfy  the  relation 

An  example  typical  of  this  behavior  is  shown  in  Figure  2. 
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Figure  2.  Divergent  velocity  for  R  =  5  m,  =  2.68  s. 


4.  Driven  Circular  Slingatron 


We  next  consider  a  class  of  shngatrons  so  fliat  the  gyration  of  the  circular  track  is 
prescribed  by  the  angular  displacement  \j/(t)  =  Xj/o  +  27t  f  t  +  Tif  t^  for  constants  f 
3nd  f .  T5^ical  numerical  solutions  to  equation  ^4)  and  equation  (5)  are  given  in 
Figure  3  with  initial  conditions  (J)q  =  0  and  =  27cf .  Two  corresponding 
calculations  of  the  angle  0  =  i(/  —  (j)  as  functions  of  time  are  shown  in  Figure  4. 
The  last  two  plote  show  that  the  slingatron  will  eventually  achieve  phase  lock  in, 
i.e.  ^  =  0  and  0  =  0,  when  the  gyration  angle  is  quadratic  in  time,  i.e. 

To  better  understand  how  the  final  lock-in  angle  is  related  to  the  given 
parameters,  we  now  rewrite  equation  (4)  in  terms  of  6  =  y/-^  and  calling 
q  =  ft  +  f  we  get 
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Velocity  (Km/s) 


r=1m  To=7r/4  ^=0.005 


—  P=0  a=l/25 

— P=l/250  a=l/25 

—  P=l/250  0=1/50 


Figure  3.  Velocity  vs.  time  P  =  f  /P,  a  =  r/R. 


r=1m  ^0=^/4  u^O.005 


(10) 


Q" +  4n\iqQ’  +  2na{^2'Kq^  — 

2Tia{^2Tiq^  +f  )co50  =  2'K{Ti\x.q^  +f  ). 

Focusing  our  attention  on  large  values  of  t ,  consider  the  following  as5anptotic 
expansion  for  0 


Substituting  this  back  into  equation  (10)  and  keeping  only  the  largest  order  term 
gives  us 


4;r ^  [a  sin(^o ) -  a//  cos  )-//]=  0  => 


However,  Tidman  [1]  has  shown  that  6  values  near  to  n/l  lead  to  instabilities, 
which  tells  us  that  the  correct  asymptotic  must  have  the  following  form 


lim—  — >  2«71  +  tan  ' 
2 


n  = 


(11) 


The  limiting  value  plotted  in  Figure  4  corresponds  to  equation  (11)  with  n  =  0 . 
Subtracting  this  asymptotic  value  of  0  so  that  y-^O  +  O^  transforms 
equation  (10)  (after  some  algebra)  to  the  following  expression: 

y'-p  f^y'^  +  4'K\xf  qy' + 2Ti{^2%Xq^  +  fp)5iny  + 

2n(2n^  q^  -f  ^)cosy  =  27t(p  ti  q^  +f  ),  (12) 

X  =  -y/(a'^-7)p-^+a^,q  =  ft  +  f. 


Since  «  //  +  /i/a  is  small  for  realistic  values  of  fj.  and  a ,  we  now  make  the 
small  angle  approximation,  y  « 1,  reducing  equation  (12)  to 


fV''  +  4;rf;£/^/  +  2;/r(2;;r2.^^  +f//)y  =  2;r  (fA  +  l) 
y(^)  =  Wo-^o  andy(f)=0. 


There  are  many  ways  to  write  solutions  to  equation  (13),  but  we  use  tiie 
following  procedure,  which  seemed  expedient  for  both  numerical  and  analytic 
examinations.  First,  make  the  substitutions 


6 


(14) 


P  =  2Tq^lf 

z  =  /exp((;r//  +  r)9^/f  ) 

into  the  left-hand  side  of  equation  (13),  which  gives  us. 


pz  + 


\ 


1 

i:-P 
2  y 


,  pU  X  +T-n:  p'-r 

z  +-^ - - — czi - z  = 

4r' 


z  =  0 . 


Next,  choose  F  =  ±i^->J  X  —  p^  to  force  the  coefficient  of  z  to  be  independent 
of  p  so  that  the  left-hand  side  of  equation  (13)  is  now  transformed  into  the 
complex  plane  such  that 


pz  +  --P 


■z/4  =  0. 


(15) 


The  solution  to  equation  (15)  is  a  hnear  combination  of  the  confluent 
hypergeometric  function  [3],  ,F, ,  so  the  solution  of  the  left-hand  side  of 
equation  (13),  Yh  /  is 


rN(p)= 

xexp 


Cl  i  pj  +  C2V:^,F,[^|,|,p^ 

i^X-p^  -p);TqYf],  X-p^  >0, 


(16) 


for  complex  constants  Cl,  C2 .  Inspecting  equation  (13)  reveals  that  the  complex 
conjugate  of  equations  (14)  will  cause  equation  (15)  to  be  replaced  by  its 
conjugate.  Thus,  the  conjugate  of  equation  (16)  will  also  produce  equally  valid 
solutions,  Yh  ,  to  equation  (13).  The  required  real  values  of  y  are  now  foimd  by 
taking  Hnear  combinations  (average  values)  of  Yh  Yh  • 

These  restdts  show  that  the  ampHtude  of  y  and  therefore  0  decrease  as  the 
factors 


(17) 


provided  that  the  functions  jFj  for  the  parameters  used  here  do  not  diverge. 
This  is  easily  shown  by  the  use  of  the  integral  representation  [3]  of  ,  F, , 

-Jjt 


,/;(1/4,1/2.a^)  = 


r(V4)r(i/4) 


(18) 


Now  taking  the  absolute  value  of  each  side  of  equation  (18)  and  remembering 
that  p  is  purely  imaginary,  we  find 
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,p;(i/4,i/2,^)|  = 


< 


r(i/4)r(i/4)| 


r(i/4)^ 


(19) 


=u 


showing  that  ,F,(l/4,l/2,p)  is  bound;  a  similar  argument  shows  that 
|,f-,(l/4.3/2,p)|  is  also  bound.  Thus,  the  homogeneous  solutions,  given  by 
equation  (16),  decay  according  to  equation  (17).  Obtaining  solutions  for  title 
limiting  case,  in  which  f  — >  0 ,  is  easily  preformed  directly  from  the  first  of 
equations  (13).  For  this  case,  the  solution  0  has  the  form 


0=  Cl  cos[27ii^a- y  f]+C2  sin\lnf-la- ^  /]  expi-lni  fit). 


(20) 


which  shows  that  the  decaying  factors  of  equations  (17)  are  replaced  with 
exp(-27tf  pt). 


The  nonhomogeneous  solution  of  equation  (13)  can  be  shown  to  asymptotically 
approach  zero  for  large  t.  Therefore,  the  asymptotic  solution  foimd  in 
equation  (11)  continues  to  show  0  «  p  +  p/a  for  small  values  of  p  and  a . 

A  sample  from  Figure  4  that  exhibits  both  the  as5nnptotic  and  decaying 
characteristics  of  a  solution  to  equation  (13)  is  given  in  Figure  5.  The  velocity 
corresponding  to  the  asymptotic  solution  of  equation  (13)  can  be  obtained  from 
equation  (5)  by  letting  ^\.  Thus,  one  will  find 


■n/  + 1 


a 


(21) 


and  this  shows  that  velocity  is  as)nmptotically  proportional  to  R  for  our  choice 
of  Iff .  Choosing  f  =  0 ,  one  will  find  diat  the  asymptotic  solution  to  equation  (10) 
wiU  now  contain  both  even  and  odd  powers  of  t ,  but  tibe  0-th  order  term  will 
still  have  the  solution  given  by  equations  (11). 

Equation  (6)  and  the  loiigtime  Hmit foimd  in  equations  (11)  tell  us  that  b  ^ 0, 
which  implies  C  =  ^  •  Thus,  the  sled  will  continue  to  gain  velocity 

for  a  long  period  of  time  whenever  the  gyration  speed  varies  quadradically  with 
t .  An  example  of  this  is  given  in  Figure  6,  which  shows  that  both  the  numerical 
and  the  analytical  results  are  in  good  agreement,  again  providing  compelling 
evidence  for  the  decaying  factors  given  in  equation  (17)  as  well  as  the  value  as 
the  final  phase-lock  angle. 

So  far,  we  have  assumed  that  the  friction  coefficient  is  very  small  but  finite 
0  <  fi  « a  <  \ .  Equations  (17)  and  (20)  show  that  friction  plays  a  very 
important  role  in  that  it  is  responsible  for  damping  out  the  oscillations  of  the 
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- Calculation  Eq.  (4) 


mass  sled  as  it  progresses  along  the  circular  track.  Setting  //  =  0  in 
equations  (11)  and  (21)  shows  that  the  long  time  limits  for  ^andV  become 
2n;rand  {pc  +  X^y/,  respectively.  Figure  7  displays  the  influence  of  friction 
where  one  can  see  that  the  velocity  no  longer  damps  to  the  asymptotic  limits 
V  =  {a  +  \)  \jf ,  when  the  force  attributed  to  friction  is  assumed  to  be  zero. 


rsim  a=0.20  %=7c/4  |i=0.0050 


- 11=0 

- jj.=0.0050 

Eg.  (21) 


Time  (s) 

Figure  7 .  Comparison  of  velocities  with  and  without  damping  p  0,  p  =  0. 


5.  Conclusions 


The  analysis  presented  here  gives  an  explicit  expression  for  the  rate  at  which  an 
accelerating  mass  becomes  phase  locked  while  it  traverses  a  circular  slingatron 
having  a  gyration  phase  angle  y/  varying  as  a  quadratic  with  time.  We  have  also 
shown  that  the  final  phase-lock  angle  has  the  simple  expression  0  — >  p  |i/a  for 
the  usual  small  values  of  p  and  a.  In  order  to  keep  0  constant,  we  must 
gradually  increase  \|/  so  that  it  diverges  at  t  =  {h%y\  as  does  the  phase-locked 
([).  This  is  because  as  the  gyration  speed  )  increases,  so  must  the 

mass  gain  speed  -27^ asin{$)  (see  equation  [5]  for  0  =  constant)  per  turn. 
However,  we  choose,  y/  =  y/Q+2;i:ft  +  7rft^,  which  causes  the  gyration  speed  to 
become  large  at  large  times.  Thus,  we  have  shown  that  the  only  way  the  mass 
can  stay  phase  locked  with  this  increasing  gyration  speed  is  to  assume  smaller 
values  of  0  so  that  it  does  not  gain  too  much  velocity  per  turn.  For  our  choice  of 
,  we  have  0  — >  p  +  p/a  «  1 ,  as  time  becomes  large,  thus  keeping  the 
accelerating  mass  phase  locked  with  gyration. 
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